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The ground-state energy of a quantum impurity model can be calculated using the numerical renormaliza- 
tion group (NRG) with a modified discretization scheme with sufficient accuracy to reliably extract physical 
information about the system. The approach is applied to study binding of magnetic adsorbates modeled by the 
Anderson-Newns model for chemisorption on metal surfaces. The correlation energy is largest in the valence- 
fluctuation regime; in the strong-coupling (Kondo) regime the Kondo-singlet formation energy is found to be 
only a minor contribution. As an application of the method to more difficult surface-science problems, we study 
the binding energy of a magnetic atom adsorbed near a step edge on a surface with strongly modulated surface- 
state electron density. The zero-temperature magnetic susceptibility is determined from the field-dependence of 
the binding energy, thereby providing an independent result for the Kondo temperature Tk, which agrees very 
well with the Tk extracted from a thermodynamic calculation. 

PACS numbers: 72.15.Qm, 73.20.Hb, 68.43.-h, 05.10.Cc 



The magnetism of nanoscopic objects supported on sur- 
faces is of great current interest due to possible applications 
in ultra-dense data storage. The magnetic properties of ad- 
sorbates can now be studied on the single-atom level using 
scanning tunneling microscopes (STM) [IJ. Adsorbed atoms 
attach to metal surfaces by forming strong (covalent) bonds 
in a process named chemisorption |0,[3t]. The chemisorption 
controls the valence (and thus the magnetic moment) of mag- 
netic adsorbates, it can lead to adsorbate-induced restructur- 
ing of surfaces, it affects superlattice growth, chemical reac- 
tions (catalysis) and other surface phenomena fstl- Using an 
STM, adsorbed atoms may be manipulated to form artificial 
nanostructures |4]. For successful manipulation of atomic- 
scale objects it is crucial to understand the binding properties 
of adsorbates, i.e. to know the potential-energy surface as a 
function of the position of the adsorbate 101 ■ 

A highly simplified model for studying the chemisorption is 
the Anderson-Newns model jslS]: = ^band + ^^imp + He 
with 

-H"band = ^ efecj-g-Cfeo-, 

k,<Te{hl} 
Hunp= ^ en^ + Un^ni, 

Hhand describes the continuum of conduction-band electrons 
with dispersion ek, i?imp corresponds to an adsorbate level 
with energy e and electron-electron repulsion U (n^ — d\d„ 
is the level occupancy operator), while iJhyb defines the hy- 
bridization which can be fully characterized by the function 
r(cj) — |Va;P(5(lj — efe). The adsorbate binding energy 
AE is defined as the difference between the ground state 
energy of the system described by the full Hamiltonian H 
and the ground state energy of the decoupled system with 
ffhyb = (i.e. the limit where the atom is far away from the 



surface). While the Anderson-Newns model was originally in- 
tended to describe binding of hydrogen and alkali atoms, some 
properties of magnetic adsorbates can also be studied within 
a single-orbital approximation [7]. General binding proper- 
ties can be determined qualitatively correctly using the unre- 
stricted Hartree-Fock method [6], while the contributions due 
to correlations can be calculated variationally A method 
which could very accurately solve the problem in full gener- 
ality for arbitrary energy-dependent r(a;) and for arbitrarily 
large interaction strength U has been, however, lacking. In 
this work, it is shown that the binding energy can be calculated 
with an excellent accurac^sing the numerical renormaliza- 
tion group (NRG) Laiail- 

The NRG consists of a logarithmic discretization of i/band 
into intervals [A"(^+^) : A^^] with A > 1, followed by a 
mapping to an effective one-dimensional tight-binding Hamil- 
tonian with exponentially decreasing hopping constants cx 
A~'/^, and an iterative diagonalisation where one further site 
is taken into account at each step. At each iteration i, the cal- 
culated excitation spectrum is shifted by subtracting the low- 
est eigenvalue Et from all others. The series 



E^KG — X El (2) 

is the ground-state energy of the effective Hamiltonian. To 
improve the results, several independent NRG calculations 
are performed for interleaved discretization meshes shifted by 
h}^^' with z S (0 : 1] and the final result is obtained as an av- 
erage over all z lfl2i [isll . To the best knowledge of the author, 
the quantity £^nrg has never been used to extract physical in- 
formation about the system, presumably due to poor conver- 
gence properties and systematic eiTors of the conventional dis- 
cretization scheme. These deficiencies of NRG were recently 
surmounted by a different discretization approach lfl4l [Tsll 
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which consists of solving the differential equation 

<i(x+l) - 



dx 



(3) 



vnx)] ' 

where the function E{x) with x = j + z yields the discretiza- 
tion coefficients for each interval j and each parameter z; the 
function e{x) defines the discretization grid [15.1 . 
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Figure 1 : (Color online) A-dependence of the calculated binding en- 
ergy of a non-interacting adsorbate. Exact binding energy 
is subtracted from the numerical results Ai5NRG(A) (circles). Full 
line is a linear fit to results in the interval A € [1.6 : 2]. The er- 
ror of the extrapolated A — > 1 value is 3.3 x 10"". The standard 
deviation ctab characterizes the spread of the results for different 
parameters z. An example of Ai5NRG(2) for A = 12 is shown in 
the inset. — 32 different values of z were used, while the pa- 
rameter iScutoff = lOuiN defines the truncation cutoff in the NRG 
iteration fT4ll . 

We first consider the binding of a non-interacting adsorbate 
with J7 = 0. In this case, the binding energy can be calcu- 
lated numerically to arbitrary precision by a simple quadrature 
(Eq. (39) in Ref. 6). For simplicity, we first consider a constant 
hybridization function: r(a>) = F for w € [—1 : 1] and zero 
otherwise. By comparing the NRG results with the exact value 
for a range of discretization parameters A, Fig.[T] we find that 
the binding energies are calculated with high accuracy even 
at A = 12; for A = 2 the error is 2 x 10"'''. If bare model 
parameters (bandwidth, e, U) are of the order of the eV, this 
magnitude of the error implies that it is possible to determine 
the binding energy with /leV accuracy. The spread of the re- 
sults Ai?NRG(2) for different values of z, as measured by the 
standard deviation (tab in Fig. [T] is not an indication of the 
error committed but rather contains physically relevant infor- 
mation about the effects of the hybridization. The z-averaging 
is thus an essential element of the binding energy calculation 
and not merely an ad-hoc procedure to accelerate the conver- 
gence. 

At large A, the error can be decreased somewhat by increas- 
ing Nz, but the improvement is minor. A systematic improve- 
ment by one order of the magnitude can, however, be obtained 
by interpolation between the data points, followed by an inte- 
gration over z on the interval [0 : 1]. The error is thereby 



reduced to 3 x 10"'' even at A = 12 with no additional cal- 
culations. (There is actually no need to use a uniform mesh 
of parameters z; it is more economical to choose the z-values 
as the quadrature nodes.) Using conventional discretization 
schemes, the eiTors are larger by orders of magnitude and 
even the extrapolated A 1 value disagrees with the ex- 
act result by 3 X 10"*; this corresponds to an eiTor of order 
meV, which is barely acceptable especially when small ef- 
fects are considered, for example in possible applications to 
long-range adsorbate-adsorbate interactions |16]. The use of 
the improved discretization scheme from Ref. 14 is thus cru- 
cial and, furthermore, the possibility of obtaining reasonably 
accurate results even at large A implies that calculations can 
be performed very efficiently. 
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Figure 2: (Color online) Binding energy AE (right vertical axis) 
and numerical error ABnrg — AiJoxact (left vertical axis) of a non- 
interacting adsorbate as a function of the hybridization F. The pro- 
portionality coefficient A_B/r = —0.70 is extracted in the interval 
F G [0 : 0.01]. For reference, the inset shows AE/T as a function 
of e. 



For large hybridization F, the adsorbate perturbs the con- 
ductance band more strongly. In NRG calculations, this is 
of particular concern since a finite representation of the band 
is used, thus finite-size effects are expected to become size- 
able. We find, however, that at A = 2 the error is bounded by 
1.7 X 10"^ for all F in the interval [0 : 0.4], Fig.|2] The bind- 
ing energy AE is linear in F to a good approximation and the 
coefficient of proportionality increases in absolute value as e 
approaches the Fermi level (see the inset in Fig.|2l), where the 
hybridization is more effective in binding the adsorbate. The 
adsorbate tends to form a bond with the substrate by sharing 
an electron with the conductance-band states. This process is 
more efficient when states in the vicinity of the Fermi level are 
involved, since their occupancy can be inexpensively changed 
by the hybridization. This is similar to bond formation in two- 
atom molecules, where the binding energy is largest when the 
atomic levels are aligned. 

We now study the full Anderson-Newns model with fi- 
nite interaction U and make comparison with the mean-field 
results obtained using the unrestricted Hartree-Fock (HF) 
method (which neglects correlation effects, see also Ref. 8). 
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Figure 3: (Color online) a) Binding energy, b) correlation energy 
Ec — AE — AEhf, Kondo temperature Tk (right vertical axis), 
and c) charge fluctuations of the single-impurity Anderson model as 
a function of the electron-electron interaction U. The Kondo temper- 
ature is extracted from the thermodynamic properties of the model 
according to Wilson's prescription fcaT/f Ximp(fcsrif )/(flMs)^ ~ 
0.07, where Ximp is the impurity magnetic susceptibility [3]. 



The binding energy reaches its highest absolute value for 
e + U ^ when the single-particle level for an additional 
electron crosses the Fermi level, Fig. [3^. This behavior is 
similar to that of the non-interacting model: the binding en- 
ergy is large when the charge fluctuates strongly. Both NRG 
and Hartree-Fock give the same qualitative features, but it is 
found that HF underestimates binding. The additional bind- 
ing energy can be defined as the "correlation energy": Ec = 
AE — AEhf- The correlation energy is largest in the valence 
fluctuation regime for e + U « F, see Fig. At this point 
the local moment begins to form (see the decreasing charge 
fluctuations Sn^ in Fig. ^ for increasing U) and the energy 
scale of magnetic correlations (the nascent Kondo regime) is 
the highest (see the Kondo temperature Tk in Fig. [3J5). The 
"Kondo-singlet formation energy" of the order of Tk does 
not account for the totality of the correlation energy: it is only 
a small fraction, in particular in the laige-U limit where the 
Kondo temperature is strongly suppressed. The most impor- 
tant contribution to the correlation energy thus stems from lo- 
cal charge correlations, rather than from extended Kondo cor- 
relations. 

The energy gain due to Kondo correlations is lost in a strong 
magnetic field, see Fig. |4] The quadratic reduction for low 
fields {gfiBB <^ ksTK) is expected due to finite spin sus- 
ceptibility at zero temperature in the strong-coupling regime 
f9]. From the prefactor we can extract the zero-temperature 
magnetic susceptibility 
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Figure 4: (Color online) Binding energy and spin polarization of a 
magnetic adsorbate in an external magnetic field (expressed in re- 
duced units of a; = gfisB /kBTK)- The inset shows the spin- 
resolved impurity spectral function in a strong field. 



where W « 1.29026 is the Wilson number H [l3]. From 
X{T ^ 0) we then obtain a value T'j^ = 1.43 x 10"^ 
for the Kondo temperature, which differs from the value of 
Tk = 1.46 X 10^* determined in a thermodynamic calcu- 
lation of magnetic susceptibility by < 3%. Considering that 
the values are obtained using entirely independent procedures, 
their close agreement is an exceptional confirmation of the 
method. The remaining small discrepancy stems mostly from 
the error associated with obtaining the coefficient of the 
contribution to the total energy in the Umit of very small mag- 
netic fields. For large fields {gfiBB ^> ksTK), the Zeeman 
effect takes over and the variation of AE is approximately 
linear. 
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Figure 5: (Color online) A-dependence of the binding energy of 
a magnetic adsorbate hybridized to a band with energy-dependent 
density of states p{uj). The hybridization function is T{lu) = 
Fopi^j) with a) sharp step function p{io) = 1 + 9{uj — ujo) 
with UJO = —0.1, b) rounded step-function p{lu) = 1 + (1/2 + 
(I/tt) tan~^ [''"('^ — '■^o)/A]), where A — 0.001, and c) oscillatory 
p{uj) = l-l-(l/2) cos[(9/2)7r(l + Lj)]. TheA = 1.6 results are used 
as reference values and subtracted from Ai5(A). 



4 



Albeit constant hybridization is a convenient simplifica- 
tion, in realistic problems r(u;) is strongly energy depen- 
dent. Three forms are considered here: sharp and rounded 
step functions, and an oscillatory function. The convergence 
with A depends significantly on the form, see Fig. |5] while 
the error remains approximately constant at « 10^® for the 
rounded step function, it increases significantly for the sharp- 
step and oscillatory function. As expected, sharp discontinu- 
ities and variations that occur over extended energy intervals 
lead to larger errors than smooth localized changes. 
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Figure 6: (Color online) Properties of a magnetic adsorbate on a 
surface at a distance of x from a step edge modelled as a hard- 
wall potential scatterer. Fitting by oscillatory power-law functions 
A + B cos{2kFX + S)/x°' provides the phase shifts <5 (show in the 
figure) and decay constants a: the binding energy AE decays as 
1/x^''^, the occupancy as 1/x^'^^, and the Kondo temperature as 
1/x^'^^. The fitting procedure was performed with the data in the 
asymptotic region x G [40 : 80]A. The inset shows the hybridiza- 
tion function at a: = 16 A. LOO = -0.039, kp = 0.217 



is thus 

Tiiu) = Ffe + Ts9 {ij - iJo) {1 - Jo {2k{ij)x]} . (5) 

While it is by now established that for magnetic impurities on 
noble-metal surfaces Tg <C lll8Ul9|] . we will nevertheless 
take a greatly exaggerated ratio Fs/F;, = 1 to accentuate the 
effect of the energy-dependenceof T{llj). In fact, on surfaces 
with giant Friedel oscillations jlSi such ratio might be realis- 
tic. 

The adsorbate properties reflect the oscillatory features in 
F(tj), see Fig. |6] The Kondo temperature is strongly coiTe- 
lated with the variation of F at the Fermi level and it can be 
well described by a cosine function with constant phase shift 
Stk multiplied by some envelope function which is, to a good 
approximation, a power-law decay The binding en- 

ergy, however, exhibits some additional structure, in partic- 
ular for low values of x. (This is not a numerical artefact: 
the same result is obtained for other choices of NRG parame- 
ters.) Qualitatively similar features are visible in the adsorbate 
level occupancy (n), but at shifted positions x. The origin of 
these effects is thus in the details of the energy dependence of 
r{u!) over a wide energy interval (i.e. on the atomic scale of e 
and U). This is unlike the Kondo temperature, which depends 
mostly on the values of T{llj) in the narrow interval on the 
scale of itself and therefore simply follows the variation 
of F((jj = 0). It may be noted that strong binding corresponds 
to high Kondo temperature and that variations of AE and Tk 
are of the same order of magnitude, pointing to a large effect 
of magnetic correlation in this situation. 

The NRG method is a very capable tool for studying corre- 
lation effects in magnetic adsorbates on surfaces. The demon- 
strated favorable scaling of errors with A brings more realistic 
(multi-orbital) models within the reach of modern computing 
facilities. The technique for calculating ground-state energies 
is very general and it can be, for example, applied to calculate 
the response of the system (expectation values, susceptibili- 
ties) with respect to arbitrary perturbations. 



The capabilities of the method for problems with strongly 
energy-dependent hybridization are demonstrated with the ex- 
ample of a magnetic adsorbate in the vicinity of a step edge 
on a surface supporting a surface-state band. The adsorbate 
hybridizes both with the bulk states via Ff, (which will be as- 
sumed not to vary with energy) and with surface states via 
Fj. On a clean flat surface, r{uj) = F;, + Ts9{u! — ujq), 
where ujq is the onset of the surface-state band. More inter- 
esting situation occurs when the adsorbate is adsorbed near 
a step edge, where the local density of states of surface- 
state electrons is modulated by standing waves. Modelling 
the step edge as a hard-wall potential, the energy-resolved 
charge density is 6n{x,u}) oc 1 — Jo[2k{uj)x], where Jo is 
the Bessel function, k{oj) the wave-number at energy u; and x 
the distance from the step edge. Modelling the surface-state 
electrons as free electrons with effective mass m* we have 

1/2 

k{uj) = [{2m* / h) {lj — cjq)] ■ The hybridization function 
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